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fT^ ■ Abstract. Although AGN feedback through ionised winds is of great importance in models of AGN/galaxy coevolution, the 

^Nj ' mass and energy output via these winds, even in the nearby universe, is poorly understood. The issue is complicated by the 

' wide range of ionisation in the winds, which means that multiwavelength observational campaigns are required to obtain the 

complete picture. In this paper, we use a ~ 160 ks XMM-Newton RGS spectrum to get the most accurate view yet of the ionised 
' outflow (warm absorber) in NGC 7469 as seen in X-rays, finding that there is a wide range of ionisation, with log ^ in the range 

fH ^ ~ 0.5-3.5 erg cm s"', and two main velocity regimes, at 580-720 and 2300 km s"', with the highest velocity gas being the least 

^ l' ionised. The total absorbing column density in the X-rays is of order 3 x 10"' cm"^. We find that the lowest ionisation phase 

of the absorber is probably identical with one of the phases of the UV absorber discovered in previous studies. We show that 
both X-ray and UV absorbers are consistent with an origin near the base of a torus wind, where matter is being launched and 
accelerated. Calculating the mass outflow rate and kinetic luminosity of all the absorber phases, we demonstrate that the X-ray 
^ I absorbing gas carries respectively ~ 90% and 95% of the mass and kinetic energy output of the ionised outflow. 
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1 . Introduction Chandra and XMM-Newton have become available. Although 

most of these analyses have been concerned with the low- 

Ionised winds from AGN are an important and somewhat enig- ^^^^^-^^ ^^^^^^^ ^^^^^ ^^^^^^^ ^^^^^-^^^ ^^^^ 

matic factor in the feedback process controlling black hole and ^^^^ signal-to-noise in the X-rays, multiwavelength studies 

galaxy evolution. Although they have an important role in a high-velocity quasar outflows, like that in PDS 456 (O'Brien 

number of popular evolutionary schemes (e.g. Silk & Rees, 2005) have begun to appear. Most attention to date has 

1998; Fabian, 1999; Granato et al., 2004; Page et al., 2004), ^^^^^^^ ^j^^jj^^j. ^j^^ and X-ray absorption lines come 

their physical properties - such as mass outflow rate and ki- ^^^^ ^j^^ ^^^^ ^^^^ ^-^^^ ^1^^ ^^^^^^ conclusion being that in 

netic luminosity - are not observationaUy weU established. The ^^^^ ^^^^^ ^^^^ ^^^^ ^^^^^ ^^^^ ^ ^^^^^-^ 

situation is complicated by flie fact that flie winds contain gas at 2005). It has been clear for a while, however, that proba- 

a wide range of ionisation levels, and so astronomers working ^jj ^-^^^ ^^^^^^^ ^p^^^j.^! signatures in both bands, 
in a single wavelength range, on either X-ray 'warm absorbers' 

or UV 'intrinsic absorbers', will receive an incomplete picture NGC 7469 is a nearby (z=0.0164, de Vaucouleurs et al. 

of the phenomenon. 1991) Seyfert galaxy that has been the target of a number of 

Combined X-ray and UV spectroscopic studies of bright observing campaigns by X-ray and UV observatory satellites, 

nearby AGN have become increasingly popular (e.g. Mathur In recent years attention has focused on the outflow of ionised 

et al., 1994, 1995; Collinge et al., 2001; Crenshaw et al., 2003; gas from its nucleus which leaves spectral traces across the X- 

Yaqoob et al., 2003; Rraemer et al., 2003), especially in the ray and UV bands. The relationship of the different velocity 

last five years since high-resolution soft X-ray spectra from and ionisation phases of the gas seen in the two wavebands, 

and the place of the outflow within the nuclear region, have 

Send offprint requests to: A. J. Blustin been of particular interest. Our previous combined study of the 

* E-mail: ajb@mssl.ucl.ac.uk X-ray and UV properties of the warm absorber in this source 
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(using a 40 ks XMM-Newton observation in 2000 and a FUSE 
spectrum obtained a year earlier Blustin et al. 2003; Kriss et al. 
2003) indicated that the X-ray absorber had a wide range of 
ionisation, and that the principal phase of X-ray absorption 
was identical with one of the two velocity/ionisation phases 
in the UV absorber. A study of simultaneous Chandra, FUSE 
and HST STIS spectra of NGC 7469 (Scott et al., 2005), ob- 
tained in 2002 (with some extra STIS observations from 2004), 
was broadly in agreement with our earlier conclusions, but nei- 
ther the Chandra HETG spectrum of Scott et al. (2005) nor the 
XMM-Newton RGS spectrum of Blustin et al. (2003) had high 
enough signal-to-noise to satisfactorily constrain the properties 
of the X-ray absorber. 

At the end of 2004, the deepest ever X-ray spectroscopic 
observational campaign on this source was performed using the 
XMM-Newton RGS. In this paper we use the resulting high 
signal-to-noise soft X-ray spectrum, together with results from 
the Scott et al. (2005) UV spectroscopic campaign, to estimate 
the relative contributions of the X-ray and UV absorbers to the 
overall mass outflow rate and kinetic luminosity of the AGN 
wind in NGC 7469. We also investigate what the new data tell 
us about the origin of the ionised outflow in this object. 

2. Data analysis 

NGC 7469 was observed by XMM-Newton for a total of 
164 ks during two long observations in November/December 
2004. During these observations, EPIC-pn was operated in 
Small Window mode, the EPIC-MOS cameras were in Timing 
Uncompressed mode and the RGS was in Spectroscopy mode. 
The OM took long sequences of 1000 s exposures in Image 
mode, starting in both cases with a cycle through the V, U, 
UVWl, UVM2 and UVW2 filters and continuing in UVW2 
until the end of the respective observations. The full observa- 
tion details are listed in Table 1 . 

This paper deals only with the RGS, EPIC-pn and OM 
data. The RGS and EPIC-pn data were reduced under SAS 
V6.1; the RGS data were processed using rgsproc, which ex- 
tracts spectra and background- subtracts them using regions ad- 
jacent in cross-dispersion space to that containing the source. 
Since the response matrices generated by rgsproc take into ac- 
count the diff"erences in effective area between the two RGS 
units, it is possible to combine both the RGSl and RGS2 spec- 
tra and response matrices (see Page et al., 2003); both first 
and second order spectra were included. Remaining instrumen- 
tal effects, not yet included in the calibration, were removed 
with the aid of residuals to power-law fits to the pure contin- 
uum source Markarian 421. The EPIC-pn data were processed 
first using epproc, and then source and background spectra and 
lightcurves were extracted from circular extraction regions of 
40" radius. For both the RGS and EPIC-pn data, the event lists 
were screened to remove intervals of high proton background 
before spectral extraction took place. The total good time in- 
tervals for the RGS and pn were, respectively, 157 ks (a 4.4% 
data loss) and 163 ks (a 0.2% data loss). 

The OM image mode photometry was done using the 
omichain pipeline running under SAS V7. The resulting source 
Usts contain count rates for each source detected, with correc- 



tions applied for various systematic efl'ects including detector 
deadtime and coincidence loss. We also corrected for Galactic 
reddening using an E(B-V) of 0.069 obtained from the NASA 
Extragalactic Database (NED), and calculated monochromatic 
fluxes using the AB magnitude system (Oke, 1974) flux con- 
version factors provided in the merged source list. 

Errors quoted for fitted or calculated quantities are for 
= 2 (the RMS of the A^^ distribution) throughout this paper. 



3. EPIC-pn spectroscopy and UV/X-ray lightcurves 

The EPIC-pn spectra were used to help define the soft X-ray 
continuum to aid in modelling the RGS spectrum. The gen- 
eral shape of the pn spectrum is similar to that observed in 
previous XMM-Newton and Chandra observations of this ob- 
ject (Blustin et al., 2003; Scott et al., 2005); there is a high- 
energy power-law with a soft excess. Firstly, the Ughtcurves 
were examined for signs of spectral variability that could affect 
the apparent form of the spectrum. Fig. 1 shows the 2120 A 
(OM UVW2), 0.2-2 and 2-10 keV pn Ughtcurves alongside 
the 2-10 keV/0. 2-2 keV hardness ratio. The spectrum hardens 
over the course of the observation; a comparison of pn spec- 
tra from the first and second halves of the observation (Fig. 2) 
reveals that although the total flux in the soft excess decreases 
between the two observations, the basic shape of the soft ex- 
cess remained the same. Detailed analysis of the variabiUty is 
presented elsewhere (Blustin et al., 2006), but for our purposes 
here, the soft X-ray continuum shape does not seem to be seri- 
ously afl'ected by spectral variability. 

In all X-ray spectral fits reported in this paper, we have in- 
corporated neutral absorption due to our Galaxy at a column of 
Nh = 4.82 X 10^'^ cm-2 (Elvis et al., 1989). Fitting a power-law 
to the 3-10 keV range, excluding 5.5-7 keV, we find a slope F 
of 1.81+0.01. In the 2000 observations, the spectral index was 
1.75 (Blustin et al., 2002), and in 2002 it had increased to 1.79 
(Scott et al., 2005). There is a Fe Ka line, weU-fitted with a 
Gaussian profile, with a FWHM of 6900+1400 km s consis- 
tent with that observed with the Chandra HETGS two years ear- 
lier (Scott et al., 2005) and with the first XMM-Newton obser- 
vation four years earlier (Blustin et al., 2003). A second weaker 
emission line, at a rest-frame energy of 6.97+0.03 keV, is most 
consistent with Fe xxvi. After fixing the parameters of the hard 
power-law and the Fe Ka and Fe xxvi lines, we added a black- 
body component to represent the soft excess. The 0.3-10 keV 
flux was 5.7 X 10"^^ erg s~^; in the previous XMM-Newton ob- 
servation in 2000 the flux was 4.4 x 10"" erg s"' in this range. 
The resulting 0.3-10 keV fit is very poor, with a reduced of 
14.35 (622 degrees of freedom). Fig. 3 shows the model super- 
imposed on the pn data, along with the residuals to the fit. The 
main cause of the bad fit is the inability of the model to repli- 
cate the form or spectral complexity of the soft excess, which 
clearly includes spectral features unresolved by the EPIC-pn. 
The next section describes detailed spectroscopy of these fea- 
tures with die RGS. 
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Table 1. Observation Identifiers, start and end dates and times, instrument modes and total exposure times. 



Obs-ID Observation start (UT) 


Observation end (UT) 


Instrument 


Mode 


Exposure (s) 


0207090101 2004-ll-30T20:48:52 


2004-12-01T20:43:52 


EPIC-MOS 1 


Timing Uncompressed (medium filter) 


84 502 






EPIC-MOS 2 


Timing Uncompressed (medium filter) 


84 507 






EPIC-PN 


Small Window (medium filter) 


84 564 






RGS 1 


Spectroscopy 


84 985 






RGS 2 


Spectroscopy 


84 980 






OM 


Image (V filter) 


I 000 






CM 


Image (U filter) 


1 000 






OM 


Image (UVWl filter) 


1 000 






OM 


Image (UVM2 filter) 


1 229 






OM 


Image (UVW2 filter) 


57 860 



0207090201 2004-12-03T01:04:39 2004-12-03T23:21:19 



EPIC-MOS 1 


Timing Uncompressed (medium filter) 


78 602 


EPIC-MOS 2 


Timing Uncompressed (medium filter) 


78 607 


EPIC-PN 


Small Window (medium filter) 


78 664 


RGS 1 


Spectroscopy 


79 085 


RGS 2 


Spectroscopy 


79 080 


OM 


Image (V filter) 


1 000 


OM 


Image (U filter) 


I 000 


OM 


Image (UVWl filter) 


I 000 


OM 


Image (UVM2 filter) 


1 220 


OM 


Image (UVW2 filter) 


53 000 





^ 2120 A : 


!i 1 1 1 1 1 1 1 1 1 1 1 1 


0.2-2 keV - 

*\ : 
1 1 1 1 1 1 1 1 I'l 1 1 1 1 1 


-V 


'' 2-10 keV - 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 I" 


:i 1 1 1 1 1 1 1 1 1 1 1 1 


(2-10 keV)/(0.2-2 keV) ~ 


5x10* 10^ 


1.5x10^ 2x10^ 2.5x10^ 



Time (s) 

Fig. 1. 2120 A (OM UVW2 filter), 0.2-2 keV and 2-10 keV band pn 
lightcurves, with hardness ratio (2-10 keV/0.2-2 keV). 



4. RGS spectroscopy 

4.1. Line identification 

The RGS spectrum contains a number of narrow absorption and 
emission line features. To assess the significance of these, we 
used a routine which fits a Gaussian profile to spectral features 
deviating from a smoothed continuum. We can use the resulting 
plot of Ax^ of the Gaussian fit against wavelength (Fig. 4) to 
estimate the statistical significance of the features; Table 2 lists 
the strongest absorption and emission features, significant at 
> 99.99% confidence (here a \Ax^\ of 16 is equivalent to 4cr 
significance for one interesting parameter). 

We identify the most prominent feature, significant at ~ 
16cr, with an absorption line blend of Fe xix and Ne ix Hear at a 
blueshift of 690 ±160 km s ' (with the rest wavelength appro- 




0.5 



1 Z 
Channel energy (keV) 



Fig. 2. EPIC-pn spectra from the whole observation (black dataset; 
750 counts per bin), the first (top dataset; 375 counts per bin) and 
second (bottom dataset; 370 counts per bin) halves of the observation 
with best fit power-law to the ranges 3-5.5 keV and 7-10 keV super- 
imposed (continuous line). The bottom panel shows the ratio of these 
datasets to the model. 



priately weighted for the relative abundances of iron and neon 
as listed by Anders & Grevesse (1989)). Several of the other 
lines correspond to transitions of H-like and He-like oxygen, 
nitrogen and carbon, and various states of sulphur, silicon and 
L-shell iron. The strongest lines (Table 2), for which we can 
obtain the most reliable wavelength measurements, are mostly 
blueshifted by a few hundred km s The blueshift of Ovin 
Lya is much higher than this, probably as a result of the bias 
introduced by the presence of the emission part of the P-Cygni 
profile. There is also evidence for a higher velocity component, 
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Fig. 4. Plot of Ax^ associated with fitting a Gaussian profile to features deviating from a smoothed continuum across the RGS spectrum of 
NGC 7469, in the rest frame of the AGN. Here a positive deviation is an emission feature, and a negative deviation implies a feature seen in 
absorption. Ax^ is equivalent in this case to the square of the number of sigma significance (here a Ax^ of 16 is equivalent to 4(t significance 
for one interesting parameter; this significance level is marked in each plot with dotted Unes). Possible identifications of features are labelled. 



over 2000 km s , in the O vii Hea and C vi Lya lines (see 
Fig 10). 

Four of the lines listed in Table 2 were also detected by 
Scott et al. (2005): Ovn f and Oviii Lya emission and Neix 
Hea and Fexx /112.832 absorption. The properties of these 
lines are broadly consistent between our analysis and theirs. 
The principal difference is that our FWHMs are less well- 
constrained due to the lower wavelength resolution of the RGS. 
The Oviii Lya emission line was found to have a velocity 
shift consistent with zero by Scott et al. (2005), whereas it is 
blueshifted by IBTO^J^jg^ km s~^ in our spectrum. The apparent 
velocity shift of this line is strongly influenced by the presence 
of O vm Lya absorption, so a change in the depth or outflow 
speed of the absorber could account for the discrepancy. The 
velocity shift of the O vii f line, which is less affected by super- 
imposed absorption, is more similar between the two analyses 
(-550 ± 160 in this analysis versus -200 ± 90 km s"\ which 



are probably consistent within the wavelength uncertainty of 
the RGS), so it seems most probable that the changing apparent 
velocity shift of O vm Lya is due to changing absorption prop- 
erties. Our Fexix/Neix Hea absorption line has an equivalent 
width roughly equal to the combined equivalent widths of the 
two Ne IX Hea components observed by Scott et al. (2005). Of 
these two components, one is blueshifted by 670 + 180 km s"' 
and one is redshifted by 860 + 100 km s~^; the velocity shift 
and FWHM of the blueshifted component are consistent with 
those for our Fe xix/Ne ix Hea line. 

4.2. lon-by-ion spectral modelling 

Taking the ions already identified as a starting point, we 
then used the slab model in SPEX 2.00.11 (Kaastra et al., 
1996, http://www.sron.nl/divisions/hea/spex/) to fit the absorb- 
ing columns of individual ions. We model the continuum un- 
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Table 2. Parameters of the strongest narrow absorption and emission lines in the RGS spectrum of NGC 7469 (significant at > 99.99% 
confidence in Fig. 4). Columns 1 to 6: /!„.„. „,eas, the measured wavelength in the (cosmological) rest frame in A (Oic<,/ is measured in the rest 
frame of our Galaxy); the line identification; Aresi, tab, the theoretical rest frame wavelength of the line; EW, the equivalent width of the line in 
mA (negative values correspond to emission lines); FWHM, the line FWHM in km s"' ; v,;,;^,, the redshift or blueshift of the line in km s"*. 



A-resr. iiieas 


Line 


A-resI, lab 


EW 


FWHM 


^ shift 


94. 7^+" '« 
23.51 ± 0.02 
22.06 ± 0.01 


N VII Lya 

OlGfl/ 

Ovnf 


24.781 
23.521 
22.101 


25^1 
29!^ 
-67 ± 15 


n+btiU 
|-|+E140 

QOQ+1070 


-570!f^j; 
-180!^^« 
-550 ± 160 


21.72 ± 0.03 


O vn i-l-r" 




-68!i^ 


2710+12™ 

^' ^ -790 




18.88 ± 0.02 


O vm Lyor 


18.969 


24 ±7 


1240 ± 910 


.J37Q+320 


17.55 ± 0.03 
15.97 ± 0.01 


Fe XVIII 
O VIII Lyj8 


17.591 
16.006 


-9 ±5 

30+* 


9ri+1440 

960 + 770 


.770+340 
' '"-460 

-610 ± 230 


14.17 ± 0.01 


Fe XVIII 


14.210 




n+2400 
"-0 


-760+'™ 

"-'"-210 


13.431 ± 0.007 


Fe xix/Ne ix Kea 


13.462'' 


61 ±6 


1710 ±330 


-690 + 160 


12.83 ± 0.01 


Fexx 


12.832 


35 ± 8 


1680!«o^ 


-160 ±310 



" No laboratory wavelength or blueshift are listed for this feature, since it is a blend of lines whose apparent average wavelength is heavily 

dependent upon the physical conditions and geometry of the emitting region. 
^ Appropriately weighted by the relative abundances of iron and neon as listed by Anders & Grevesse (1989). 




d 0.5 1 a 5 10 

Energy (keV) 




0.5 1 2 5 10 

Energy (keV) 



Fig. 3. Top: EPIC-pn spectrum (black) with fitted model (red) com- 
prising a Galactic-absorbed power-law, Fe Ka and Fe xxvi lines, and 
a blackbody soft excess. Bottom: the residuals of the fit above; the 
model cannot explain the spectral complexity in the soft excess. 



derlying the RGS spectrum with a hard power-law and a soft 
black-body component, as seen in the pn spectrum; the param- 
eters of this model are given in Table 3. We consider this to be 
a convenient parameterisation rather than a physical model; the 
nature of the soft excess in AGN is still under debate (see e.g. 
Brocksopp et al., 2006; Crummy et al., 2006; Gierlinski & 



Done, 2006). The velocity shift and turbulent velocity asso- 
ciated with the absorption lines were allowed to vary during 
the fit, and were assumed to be the same for all states. Since 
the deepest absorption lines for some of the most important 
species form part of P-Cygni profiles, it was also necessary to 
include the set of emission lines corresponding to these transi- 
tions; the parameters of these lines were fitted simultaneously 
with those of the corresponding absorption features. The pa- 
rameters of the slab best-fit model are listed in Table 3, and the 
fitted parameters of the emission lines are given in Table 4. The 
model is shown overplotted on the RGS spectrum in Fig. 5. The 
P-Cygni line profiles and very high fluxes of the narrow emis- 
sion lines may already imply an extended flow and a significant 
mass outflow rate, as also suggested by Behar et al. (2003) for 
NGC 3783. 

In Fig. 6, the absorbing columns of the ions in this model 
are plotted against the ionisation parameter of greatest abun- 
dance of that ion. The plot displays the usual trend (c.f. Blustin 
et al., 2005) of higher absorbing columns at higher levels of 
ionisation: this is due to an observational bias in that highly 
ionised gas needs to have a higher column in order to leave 
a visible trace in the spectrum. The lower-ionised states of 
iron fit the apparent Unresolved Transition Array (UTA) fea- 
ture between 16-17 A. Unusually, there is no evidence of 
the lowly ionised oxygen states that generally accompany iron 
UTAs (e.g. Blustin et al., 2002). If these UTA ions are indeed 
present - the slab fit improves by a of 82 when the UTA 
ions are added to it - then this may be an indication of a high 
relative abundance of iron in the lowest ionisation levels visi- 
ble to us. The fitting of the UTA ions is also dependent upon 
the assumed outflow velocity, as the spectral resolution does 
not allow the identification of individual ions within it. The fit- 
ted RMS velocity of the gas, 97 + 9 km s ', corresponds to 
a Doppler width of b - 137 + 13 km s ^ This is close to the 
value inferred as most probable by Scott et al. (2005), who esti- 
mated ^7 ~ 100 km s ' for the soft X-ray absorption lines on the 



6 



A. J. Blustin et al.: The mass-energy budget of the ionised outflow in NGC 7469 




basis of the widths of nine-velocity-component UV absorption 
features measured in FUSE spectra. 

4.3. Absorption measure distribution 

Fig. 6 shows that there is a broad range of ionisation present in 
the absorber, which spans about four orders of magnitude in ^. 
Following the method of Holczer et al. (2007), we express the 
total hydrogen column density Nfj along the line of sight as an 
integral over its distribution in log ^. We term this distribution 
the Absorption Measure Distribution (AMD - dNn /d(log^)). 
The column density of each ion Nion can be expressed in terms 
of Nh as: 



Nion ^ A, AMDfU\ogOd(log^) 



(1) 



where Ar is the elemental abundance with respect to hydro- 
gen, assumed to be constant throughout the absorber, and 
is the fractional ion abundance with respect to the total abun- 



® 


Ne 


•k Fe 


A 





■ S 


O 


N 


Si 


□ 


C 


• Mg 




Log fmax 

Fig. 6. Plot of the absorbing columns of ions in the slab model against 
the log ionisation parameter ^ for which the ion is most abundant. 



dance of its element, which we obtain from XSTAR (Kallman 
& McCray, 1982) using the SED described in Section 4.4. 
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Table 3. Parameters of the best-fit slab model for the RGS spectrum 

of NGC 7469: the slope T and normalisation Normpi (in 10" pho- 
tons s"' keV"') of the hard power-law; the temperature t (keV) and 
normalisation NormsB (10'* m^) of the soft blackbody; the average 
velocity shift v,;,,/, and RMS turbulent velocity v,urb of the absorber 
(both in km s"'); the log of the column densities of ionic species in 
the spectrum in cm~^. 



~ log(T) [K] 

4.5 5.0 



r 


1.81 + 0.04 


Norm PL 
t 

NomiBB 


5.9 ±0.1 

n 144+0.007 
84+30 


Vshift 
Vturb 


-690 ± 40 
97 ±9 


Cvi 
Nvi 
Nvn 


16 7 ^"-^ 
15 4 

^-'•^ -1.4 

16.5 ± 0.2 


Ovn 
Ovm 


16 3 ■""•^ 
17 96 


Neix 
Nex 
Mgxi 
Mgxn 
Six 
Si XI 


17 2 
17 3 

^ ' •-' -0.4 

16.4 
16.5 
16.8 
15.9 


Sixn 


16.7 ± 0.2 


Sx 

S XII 


16 2 -^-^ 
15.9 ± 0.3 


S XIII 
Feix 
Fex 
Fexi 

Fexii 
Fe XIII 


15 7 +"-3 
15 7 ' 
15.8 to,2 
15 5 

15.0 

15.7 ± 0.2 


Fexiv 
Fexv 
Fexvi 
Fe XVII 
Fexvin 


15 4 ^"-^ 

^-^•^ -0.6 

15.3 

15 3 ^o-' 

16 ^"^ 

iu.u Q3 

16.3 ±0.1 


Fexix 


16.4 ±0.1 


Fexx 


16.3+0.1 


Fexxi 


16.6 ± 0.2 


Fe XXII 
Fexxin 


15 5 "-"-^ 

16 5 



Table 4. Parameters of the emission lines included with the slab and 
xabs model fits to the RGS spectrum of NGC 7469: line identification, 
FWHM in km s~' , velocity shift VsMft in km s"' and Une flux F in 10~* 
photons cm~^ s"'. 



Line 



FWHM 



shift 



C VI Lyor 


n+2oOo 
"-0 


"'"-470 


30± 10 


N vn Lya 


Q+UOOO 


1 100+540 


9±7 


Ovnf 


1200+""'" 

^^""-700 


-520 ± 200 


70 ± 10 


Ovn i 


2500+1500 


-760+'*° 

' ""-720 


60 ±20 


vnr 


070+1800 
=""-970 


.110+1300 


40± 10 


Ovm Lyor 


640+630 




80!f» 


Ne X Lya 


n+1000 

"-0 


300+440 

-'""-230 


18 ±6 


Fexvmil7.56A 


r,+ 1900 
"-0 


-770 ±310 


15 ±6 


Fexvnil7.14A 


n+5400 
"-0 


1700+'""' 


22 ±5 




log(^) [erg cm s"^] 



Fig. 7. The Absorption Measure Distribution (AMD) of NGC 7469 re- 
constructed from the ionic column densities of iron listed in Table 3. 
The AMD has a minimum at 4.5 K< log T < 5 K, demonstrating that 
the absorber is divided into two main phases. The integrated column 
density up to ^ is plotted in the lower panel. The corresponding tem- 
perature scale is shown at the top of the figure. 



We find the AMD that best reproduces the whole set of Fe 
ionic column densities simultaneously. This distribution is pre- 
sented in Fig. 7; we get a total equivalent hydrogen column 
density of A^h = (3.3±0.8)xlO^' cm"^. More details on the 
AMD reconstruction method and the associated error computa- 
tions can be found in Holczer et al. (2007). In Fig. 7, there is a 
minimum at 1 .0 < log ^ < 2 corresponding to 4.5 < log T < 5 
(K), and a much higher column density in the high-ionisation 
component than in the low-ionisation component. The AMD 
minimum at 4.5 < log T < 5 (K) is seen in other AGN outflows 
as well (Holczer et al., 2007) and we speculate that it is due to a 
ubiquitous thermal instability in these winds. The abundances 
are obtained by rescaUng the AMD (Fig. 7) to reproduce the 
measured non-Fe ionic column densities. The resulting values 
are listed in Table 5 and are broadly consistent with Solar. 



4.4. Multi-phase warm absorber model 

Taking the ionisation distribution revealed by the AMD as a 
starting point, we constructed a multi-phase spectral model of 
the warm absorber using the xabs model in SPEX. Xabs is 
based on a grid of XSTAR output, and generates absorption 
features due to photoionised gas at a given ionisation parame- 
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Table 5. Abundances with respect to Solar (Anders & Grevesse, 1989) 
derived from the Absorption Measure Distribution (AMD. 

Element Abundance 



±1 "•'^-0.4 

O 1.3 ± 0.5 

Ne 3*1 
Mg l.O-l^ 

Si 11+9 

s 6:i° 

Fe 1^ 

" Fixed 

ter and column density. The ionisation parameter ^ is defined 
as 

^=^. (2) 

in which L,o„ is the 1-1000 Rydberg ionising source lumi- 
nosity (in erg s"'), n the gas density (in cm"^) and r the source 
distance in cm, so ^ has the units erg cm s"^ (Tarter et al., 1969). 

The Spectral Energy Distribution (SED) that we used in 
calculating the xabs models is based on that used by Scott et al. 
(2005) to represent the ionising continuum of NGC 7469 as 
it was in 2002, but adjusted for the broadband spectrum of 
the source as observed with XMM-Newton two years later. 
The hard power-law, above 2 keV, has F = 1.86 to match the 
slope fitted to the EPIC-pn spectrum between 2-10 keV (with 
6-7 keV ignored), and is normalised to the unabsorbed flux of 
the XMM-Newton spectrum at 2 keV. The soft excess begins 
to rise from ~ 2 keV, so we take this energy to define the high- 
energy end of a power-law joining the UV and soft X-rays, 
and, following Scott et al. (2005), we set the low-energy end 
at 50 eV. Below this energy, we use the same spectral form as 
Scott et al. (2005), with E = 2 below 2500 A and F = 1.92 
between 2500 A and 50 eV. The lowest-energy section was 
normalised to the V-band (5430 A) flux from the OM, and the 
resulting 50 eV-2 keV slope was calculated to be F - 2.78. 
The spectrum cuts ofFbelow 0.1 eV and above 100 keV. Fig. 8 
is a plot of the resulting SED over the spectral range of interest, 
with monochromatic fluxes for the OM photometry points and 
the pn at 2 keV overplotted. The SED of Scott et al. (2005) is 
also plotted for comparison. 

The AMD shows (Fig. 7) that the absorber contains two 
principal ionisation ranges: a low ionisation phase with log ^ ~ 
0.5 and a higher ionisation phase with an ionisation range of log 
^ ~ 2.5-3.5. We constructed a spectral model with three xabs 
components to reproduce this, with one component for the low 
ionisation gas and two components to represent the high ionisa- 
tion range. The value for the turbulent velocity obtained in the 
slab fit, 97 km s"', was used for all phases, and the abundances 
were fixed at Solar (Anders & Grevesse, 1989). The emission 
lines were included with parameters fixed at those previously 
fitted alongside the slab model (see Table 4). The absorbing 
columns, ionisation parameters and outflow speeds of the three 
phases obtained from the xabs fits are listed in Table 6. A plot 
of the spectrum with the model overlaid is given in Fig. 9. 




14 15 16 17 IB 

Frequency (Hz) 



Fig. 8. The NGC 7469 spectral energy distribution at the time of the 
XMM-Newton observation in late 2004 (solid line), as inferred from 
the OM and pn monochromatic fluxes (stars). The SED used by Scott 
et al. (2005) to represent data taken two years earlier is also plotted 
for comparison (dashed line). 

Table 6. Parameters of an xabs model fit to the RGS spectrum of 
NGC 7469, using three phases of gas at different levels of ionisa- 
tion. For each phase, the following parameters are listed: log ioni- 
sation parameter ^ (erg cm s"'), equivalent hydrogen colunm density 
Nb (cm"^), and average velocity shift Vj/,,/, in km s"'. 



Phase 


Parameter 


Value 


Component 1 


Log^ 


0.8:- 




N„ 


3!^ X lO" 




Vshift 


-2300 + 200 


Component 2 


Log^ 


2.73 + 0.03 






2.0 ± 0.2 X 10^1 




Vshifl 


-720 ± 50 


Component 3 


Logf 


3 56+'"'*' 

-^••^ -0.07 




Nh 


2.9+°-? X 10^' 




Vshift 


-5'80:»« 



Fixed 



We found that whilst the higher ionisation phases have fit- 
ted outflow velocities of 720 and 580 km s"\ the lowest ion- 
isation material was consistent with a higher velocity outflow 
phase at 2300 km s"^ . Fig 10 contains plots of the two velocity 
components associated with the Ovn Heor, Nvn and Cvi Lyor 
lines with the model overlaid. 

The low ionisation component 1 should give rise to absorp- 
tion by states of M-shell iron forming an Unresolved Transition 
Array. With Solar elemental abundances, our model does not 
reproduce the full depth of this absorption; this is probably 
due to uncertainties in the ionisation balance of the UTA. A 
recent re-calculation of the dielectronic recombination rates 
for Fe**+-Fe'^+ with a modern atomic physics code (Badnell 
2006) has indeed shown that the total recombination rates for 
these ions, radiative plus dielectronic, are an order of magni- 
tude higher than those currently used in XSTAR and other pho- 
toionised plasma codes. 
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Fig. 9. RGS spectrum composed of data from the whole observation with the three-phase spectral model superimposed: the model parameters 
are listed in Table 6. 



5. Discussion 

5.1. The mass-energy budget of the outflow 

When we want to quantify the mass transfer involved in AGN 
feedback in the context of galaxy evolution, how important is it 
to take the whole ionisation range of the outflow into account? 
To answer this, we need to know how much mass is coming out 
in all of the ionisation phases that we observe. We can estimate 
the mass outflow rates M of the UV and soft X-ray absorbing 
outflow phases in NGC 7469 using Eqn. 18 in Blustin et al. 
(2005): 



M 



\.23mpLionCyvQ. 



(3) 



an outflow consisting of, for example, narrow streamers of gas 
(Steenbrugge et al., 2005). 

We can obtain a simple estimate of the volume filling fac- 
tor Cy if we assume that the outflow is radiatively accelerated. 
In such an outflow, the momentum of an outflowing phase (de- 
pendent upon C,,) must be of the order of the momentum of 
the radiation that it absorbs. Pahs, plus that of the radiation that 
it scatters, Pscaii- The resulting expression for Cy is given by 
Eqn. 23 in Blustin et al. (2005): 



1.23to„L,o„v2Q 



(4) 



(note that the factor of c in the denominator of this equation as 
originally printed was a typo). We can write 



where nip is the proton mass, L,o„ is the ionising 
(1-1000 Ryd; 13.6 eV-13.6 keV) luminosity, C, is the vol- 
ume filling factor of a given ionisation phase, v is the outflow 
speed, Q. is the solid angle subtended by the outflow and ^ is 
the ionisation parameter (Eqn. 1). As in Blustin et al. (2005), 
we take Q.- 1.6, which is estimated from the observed fraction 
of type-1 AGN in the nearby universe, ~ 25% (Maiolino et al., 
1995), and the observation that the covering factor of warm 
absorber outflows in type-1 AGN is at least 50% (Reynolds, 
1997). Note however that Q. could be a lot smaller than this in 



abs — 



(5) 



in which c is the speed of light, and Labs is the luminosity ab- 
sorbed by the outflow over the whole 1-1000 Ryd range; the 
absorbed luminosity is obtained using the xabs warm absorber 
model in SPEX, taking the non-unity covering factors of the 
UV phases into account. Also 



C 



0, 



(6) 
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-4000 -3000 2000 4000 

Velocity (km s~^) 

Fig. 10. Velocity spectra of the O vii Hea, N vii and C vi Lyo- lines, 
showing the two principal velocity phases in the outflow. The dotted 
lines mark the fitted outflow speeds of the three xabs components, at 
2300, 720 and 580 km s ' for components 1, 2 and 3 respectively (see 
Table 6. 



where tt (the optical depth for Thomson scattering) is given 
by 



Tt = cTj Nh , 



(7) 



o-j is the Thomson cross-section and Nh is the absorbing col- 
umn. 

Having obtained the mass outflow rates, we can also calcu- 
late the kinetic luminosities of the outflow phases, using 



1 . , 

Lke = ^Mv"- 



(8) 



To get an accurate figure for the total mass outflow rate, we 

also need to know whether any of the UV absorption is due to 
the same phases of gas that are absorbing in the X-rays. We can 
test whether this is the case by comparing low-ionisation ion 



columns predicted by our xabs models with the columns actu- 
ally measured using transitions in the UV spectra. In Table 7, 
we list C IV, N V and H i column densities predicted by our mod- 
els alongside values obtained from the June 2004 STIS spectra 
(Scott et al., 2005), which are the measurements closest in time 
to our December 2004 X-ray spectra. We also list the ionisa- 
tion parameters and velocities of the two UV and three X-ray 
phases; the UV phase ionisation parameters were derived us- 
ing spectra taken in 2002, so we have rescaled them to take 
account of the greater ionising luminosity at the later epoch 
(2. 19 X 10^4 in 2004 versus 1 .45 x lO"*"* erg s"' in 2002). The ve- 
locities of the UV phases did not change significantly between 
the 2002 and 2004 observations; using the components of the 
Lyo- absorption line as an example, the column-weighted av- 
erage high- velocity component was -1896 ± 3 km s ' in 2002 
and -1893 + 4 km s~^ in 2004, whilst the low- velocity compo- 
nent was -560 + 10 km s'^ and -580 ± 10 km s'^ in 2002 and 
2004 respectively. 

The high-speed, low ionisation X-ray component 1 has an 
ionisation parameter consistent with that of UV component 2. 
Its outflow velocity is similar, especially taking the system- 
atic wavelength uncertainty of the RGS (at 100-200 km s"^) 
into account. The columns of Civ, N v and H i predicted to be 
present in this phase are very close to the actual values mea- 
sured in UV spectra. We therefore identify UV component 2 
with X-ray component 1. Of the other ionisation phases, the 
outflow speeds of X-ray components 2 and 3 may be consistent 
with that of UV component 1, but the ionisation levels of the 
phases, and predicted columns, are not consistent. This is the 
opposite conclusion to that reached by Blustin et al. (2003) and 
Kriss et al. (2003), namely that the high-ionisation phase was 
identical with UV component 1 , but UV component 2 lacked an 
X-ray analogue. We note, however, that our new X-ray spectra 
are vastly superior to the original observations of Blustin et al. 
(2003). 

Table 8 gives the observed parameters and calculated per- 
centage volume filUng factors, mass outflow rates, and kinetic 
luminosities for the X-ray and UV absorbing phases. Taking 
all of the X-ray and UV phases into account, but assuming that 
UV component 2 is identified with X-ray component 1 (and us- 
ing the mass outflow rate of the X-ray phase, which is higher), 
we obtain a total mass outflow rate of 0.06 Mg yr ' . Calculated 
in the same way, the total kinetic luminosity of the outflow is 
1 .2 X 10"*" erg s"' . The X-ray absorbing gas is responsible for ~ 
90% of the mass outflow rate and 95% of the kinetic luminosity 
of the outflow. 

In the case of UV component 2/X-ray component 1, the 
mass outflow rate and kinetic luminosity calculated from the 
UV properties are respectively 21% and 14% those derived 
from the X-ray properties. This is partly due to the slightly 
higher outflow velocity measured from the X-ray spectrum (re- 
sulting from the larger wavelength uncertainty of the X-ray 
measurements), but also to the order of magnitude lower col- 
umn density obtained from the U V. What we have not yet tested 
is how much mass and energy is transported by matter at much 
higher ionisation levels, containing ionisation states that pri- 
marily absorb at energies above 6 ke V. This will be easiest to do 
for high-column warm absorbers, Uke for instance NGC 3783, 
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Table 7. A comparison of some properties of the X-ray and UV absorbing outflow phases in NGC 7469: absorber phase, log ionisation 
parameter (log ^; erg cm s"'), outflow speed (v; km s"'), Crv colunm density (10'* cm"^), N v column density (10''* cm"^), Hi column density 

(10''' cm-2). 



Phase 


Log^ 


V 


Nciv 


Nnv 


Nh, 


X-ray component 1 




-2300 + 200 


1.6 


3.4 


6.2 


X-ray component 2 


2.73 ± 0.03 


-720 ± 50 


Negligible 


0.00091 


Negligible 


X-ray component 3 


3 56+°"^ 


-580+«« 


Negligible 


Negligible 


Negligible 


UV component 1 


1.6P 


-562" ± 6 


0.98 ± 0.09 


2.9 ± 0.8 


7±2 


UV component 2 


0.5P 


-1901" ± 6 


2.0 ±0.1 


2.5 + 0.2 


2.4 + 0.5 



" Converted from U to ^ using the Scott et al. (2005) SED, and rescaled to take into account the change in ionising luminosity between the 

2002 and 2004 observations. 
" As measured from 2002 FUSE and STIS data (Scott et al., 2005). 



Table 8. Physical properties for each warm absorber phase in NGC 7469: absorber phase, log ionisation parameter (log ^; erg cm s"'), log 
column density (log A'h; cm"-), outflow speed (v; km s"'), covering factor (/), percentage volume filUng factor (% Cy), mass outflow rate per 
phase (Mout; Mq yr '), log kinetic luminosity of outflow (log Lke', erg s"'). 



Phase 


Log^ 


LogNH 


V 


/ 






Log Lke 


X-ray component 1 


0.8 


19.5 


2300 


1 


0.0005 


0.002 


39.6 


X-ray component 2 


2.73 


21.3 


720 


1 


2 


0.03 


39.7 


X-ray component 3 


3.56 


21.5 


580 


1 


12 


0.02 


39.4 


UV component 1 


1.61^ 


20.0 


562 


0.53 


0.04 


0.006 


38.7 


UV component 2 


0.5P 


18.6 


1901 


0.93 


0.00007 


0.0004 


38.7 



" Converted from U to ^ using the Scott et al. (2005) SED, and rescaled to take into account the change in ionising luminosity between the 
2002 and 2004 observations. 



where there will be more absorption features visible at the high- 
est energies (see e.g. Reeves et al., 2004). 



We note here that although our estimates are based on the 
assumption of a radiatively-accelerated outflow, the nature of 
the acceleration mechanism in AGN winds remains contto- 
versial. Models of radiatively-driven winds in AGN have the 
difficulty that the gas can become so highly ionised that the 
optical depth is no longer great enough for radiative accel- 
eration to take place. This problem can be overcome where 
parts of the outflow are shielded from the soft X-rays by either 
a highly ionised 'blanket' layer, or so-called hitch-hiker gas 
which is caught up in the accelerating outflow (see e.g. Murray 
et al., 1995), a rather more complex scenario than the simpli- 
fied analysis we present here. Other models involve the gas 
being thermally driven (e.g. Krolik & Kriss, 2001) or flowing 
out from an accretion disc along (locally) open magnetic field 
fines (e.g. Bottorff et al., 2000). In some models this magnet- 
ically launched gas is then radiatively accelerated (e.g. Konigl 
& Kartje, 1994; Everett, 2005). Everett & Murray (2006) 
and Das et al. (2006) attempt to model the observed outflow- 
ing Narrow Line Region (NLR) dynamics in nearby Seyfert 
2s using, respectively, thermally and radiatively driven winds. 
Neither model could reproduce the acceleration profile of the 
outflowing clouds, and it seems likely that interaction with the 
surrounding medium has a significant efi'ect on the dynamics 
of AGN winds. 



5.2. The ionised outflow in NGC 7469: a torus wind 
being launctied? 

We can use the measured properties of the warm absorber in 
NGC 7469 to estimate the distances of the different absorb- 
ing phases from the central engine. Upper and lower limits for 
these distances can be obtained following Eqns. 27 and 32 of 
Blustin et al. (2005). The nunimum distances are obtained on 
the assumption that the outflowing gas must have reached the 
escape velocity, and the maximum distance from the criterion 
that the gas must be close enough for it to be ioiused. The 
minimum distance is likely to be the weaker constraint. We 
may be observing the gas before it has reached its final veloc- 
ity, so it could be closer to the centtal engine than we predict. 
Observations of outflowing gas in the NLRs of local Seyfert 2s 
shows that the wind accelerates up to a distance of order 100 pc 
(Ruiz et al., 2005; Das et al., 2005, 2006), so we can certainly 
say that it has effectively escaped from the gravitational poten- 
tial well of the black hole. It does not escape from the galaxy's 
potential weU, however, as these same observations show that 
the gas decelerates again by the time it reaches a distance of 
300-400 pc. This deceleration is very likely, however, to be 
at least partly the result of interactions with the surrounding 
medium. 

Table 9 fists the calculated distance ranges, and Fig. 11 
shows these distances plotted as ratios to the distance of the 
torus. Using dust reverberation, Suganuma et al. (2006) show 
that the iimer radius of the torus in NGC 7469 is in the range of 
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Table 9. Distance ranges of the ionised outflow phases from the cen- 
tral engine (R, pc). 



Phase R 

X-ray 1 0.012-1.7 

X-ray 2 0.13-1.3 

X-ray 3 0.20-0.81 

UV component 1 (2004 epoch) 0.21-6. 1 

UV component 2 (2004 epoch) 0.018-3.6 







1 1 


1 1 1 1 1 

irv 2 

T 

1 


1 1 1 i 

X-ray 2 

1 
1 


1 1 1 1 1 1 

X-ray 3 

1 " 






1 




1 
1 

J. 


1 : 

_L : 








_L 

X-ray 1 




^BLH 







1 


2 


3 


4 



Fig. 11. Locations of the X-ray (dashed lines) and UV (dotted lines) 
absorbing phases with respect to the BLR and torus. Y-axis: upper and 
lower bounds to the distances of the absorbing phases from the central 
engine, expressed as ratios to the distance of the torus. The ratio of the 
BLR to torus distances is also plotted (X-axis is arbitrary units). 



65-87 light-days from the nucleus. For our purposes, we take 
the torus distance to be the average of these values, which cor- 
responds to 0.064 pc. The location of the BLR is also plotted 
in Fig. 11, where we used the reverberation mapping distance 
of 5 light-days measured by Wandel et al. (1999). 

With the estimates we obtain here, we find that all of the 
phases are further out than the BLR, and most are further out 
than the torus. This contradicts the previous finding that UV 
component 1, which was identified with the high ionisation X- 
ray absorber, is probably at or within the BLR (Kriss et al., 
2003; Blustin et al., 2003). This inconsistency was first noticed 
by Scott et al. (2005), who also applied the distance estimates 
presented by Blustin et al. (2005), and discussed a number of 
reasons for the disagreement. The main evidence for the X- 
ray/UV component 1 absorber phases being located at or within 
the BLR was that UV component 1 has a very low covering 
factor of ~ 50%, consistent with a partial covering of both the 
broad lines and continuum, or covering only the UV continuum 
and not the broad Unes. Its relatively high ionisation parame- 
ter and column, which associated it with the main body of the 
X-ray absorber, supported this conclusion. Component 2 has 
a much higher covering factor of ~ 90%, which implies that 
it is further away from the central engine. The faster outflow 
speed of this more distant component indicates that the wind is 
accelerating. 



The new information we can now add to this, with much 
better quality X-ray data, is that the low ionisation X-ray ab- 
sorption is likely to be identified with the high velocity UV 
component 2, due to the similar outflow velocities and pre- 
dicted column densities of low-ionisation ionic states observ- 
able in the UV. The high ionisation phases of the X-ray ab- 
sorber are unUkely to be identified with UV component 1, prin- 
cipally because the ionic columns measured earlier in 2004 us- 
ing STIS are not predicted by our xabs models for these phases, 
even though the velocities of the phases are probably consistent 
(considering the 7 mA absolute wavelength uncertainty of the 
RGS, which is equivalent to -100 km s"' at 20A). 

The low covering factor of UV component 1 could be con- 
sistent with our distance estimates (Table 9, Fig. 1 1) if it com- 
prises a filamentary structure, perhaps near the base of a wind 
where there is acceleration of clumpy, inhomogeneous matter 
from the dusty torus. The very low volume filling factor we de- 
rive for this phase, at 0.04%, supports this interpretation. The 
X-ray absorber, at a similar speed to component 1, would also 
be located close to the base of the outflow, where one would 
expect gas at a range of ionisation levels to exist as matter is 
ionised and driven outwards. UV component 2 probably rep- 
resents gas which has been accelerated away from the base of 
the outflow, as it is moving faster, and its higher covering factor 
indicates greater distance. The X-ray and UV properties of the 
NGC 7469 warm absorber may, therefore, be the signature of a 
torus wind in the process of being launched. 

6. Conclusions 

We use our ~ 160 ks soft X-ray spectrum of NGC 7469, the 
highest quality soft X-ray data yet obtained for this source, 
to obtain accurate measurements of the ionisation state, ab- 
sorbing column and outflow speed of its warm absorber. We 
show that there is gas over a wide range of ionisation (log 
^ ~ 0.5-3.5), and we find evidence for two separate velocity 
regimes at 580-720 and 2300 km s"^ Contrary to expectations 
from previous observations, the properties of the high ionisa- 
tion X-ray absorber are not consistent with those of either of 
the components of the UV absorber, but the lowest ionisation 
X-ray absorber is probably identical with the highest velocity 
UV phase. We estimate the distances of the X-ray and UV ab- 
sorbing phases from the central engine, finding that they are 
most consistent with a location close to the torus or between 
the torus and BLR, and discuss a scenario where the cover- 
ing factors and velocities of these absorber phases signify the 
launch and acceleration of a wind from the torus. 

Using the observed properties of the X-ray absorbing 
phases, along with the luminosity-scaled properties of the UV 
absorber obtained from observations two years earlier (Scott 
et al., 2005), we calculate the mass outflow rates and kinetic 
luminosities of all of the soft X-ray and UV warm absorber 
phases, taking into account the identification of an X-ray phase 
with a UV phase. We estimate that ~ 90% of the mass outflow 
rate and ~ 95% of the kinetic luminosity are associated with the 
X-ray absorbing phases. A complete picture of the mass/energy 
budget of the multi-ionisation-phase warm absorber also needs 
to include the highest ionisation gas, which primarily absorbs 
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above 6 keV. The properties of this high-ionisation gas are most 
easily measured in outflows with higher absorbing columns 
than NGC 7469. 

For this nearby AGN, we have demonstrated that the mass 
and energy outputs are probably dominated by the X-ray warm 
absorber, although an appreciable fraction is only detectable 
in the UV band. It will be interesting to see how far this con- 
clusion extends to AGN with cosmologically-interesting BAL- 
type outflows, but it is clear that any observational study of 
AGN feedback via ionised winds needs to take the whole ioni- 
sation range of the outflow into account. 
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